Near linear time algorithm to detect community structures in large-scale networks 
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Community detection and analysis is an important methodology for understanding the orga- 
nization of various real-world networks and has applications in problems as diverse as consensus 
formation in social communities or the identification of functional modules in biochemical networks. 
Currently used algorithms that identify the community structures in large-scale real-world networks 
require a priori information such as the number and sizes of communities or are computationally 
expensive. In this paper we investigate a simple label propagation algorithm that uses the network 
structure alone as its guide and requires neither optimization of a pre-defined objective function nor 
prior information about the communities. In our algorithm every node is initialized with a unique 
label and at every step each node adopts the label that most of its neighbors currently have. In 
this iterative process densely connected groups of nodes form a consensus on a unique label to form 
communities. We validate the algorithm by applying it to networks whose community structures 
are known. We also demonstrate that the algorithm takes an almost linear time and hence it is 
computationally less expensive than what was possible so far. 

PACS numbers: 89.75.Fb, 89. 75. He, 87.23. Gc, 02. 10. Ox 

I. INTRODUCTION 

A wide variety of complex systems can be represented as networks. For example, the World Wide Web is a network 
of webpages interconnected by hyperlinks; social networks are represented by people as nodes and their relationships 
by edges; and biological networks are usually represented by bio-chemical molecules as nodes and the reactions between 
them by edges. Most of the research in the recent past focused on understanding the evolution and organization of 
such networks and the effect of network topology on the dynamics and behaviors of the system [J, Q, Hi El- Finding 
community structures in networks is another step towards understanding the complex systems they represent. 

A community in a network is a group of nodes that are similar to each other and dissimilar from the rest of the 
network. It is usually thought of as a group where nodes are densely inter-connected and sparsely connected to other 
parts of the network [J, [H, [6[ . There is no universally accepted definition for a community, but it is well known that 
most real-world networks display community structures. There has been a lot of effort recently in defining, detecting 
and identifying communities in real- world networks 0, 0, H, H, EH EH, EH EH EH EH- The goal of a community 
detection algorithm is to find groups of nodes of interest in a given network. For example, a community in the WWW 
network indicates a similarity among nodes in the group. Hence, if we know the information provided by a small 
number of webpages, then it can be extrapolated to other webpages in the same community. Communities in social 
networks can provide insights about common characteristics or beliefs among people that makes them different from 
other communities. In bio-molecular interaction networks, segregating nodes into functional modules can help identify 
the roles or functions of individual molecules [lo| . Further, in many large-scale real-world networks, communities can 
have distinct properties which are lost in their combined analysis [l| . 

Community detection is similar to the well studied network partitioning problems [l6|, EH EI| ■ The network parti- 
tioning problem is in general defined as the partitioning of a network into c (a fixed constant) groups of approximately 
equal sizes, minimizing the number of edges between groups . This problem is NP-hard and efficient heuristic methods 
have been developed over years to solve the problem [TH, EH EM fl9l-l20j] . Much of this work is motivated by engineering 
applications including very large scale integrated (VLSI) circuit layout design and mapping of parallel computations. 
Thompson [2l| showed that one of the important factors affecting the minimum layout area of a given circuit in a 
chip is its bisection width. Also, to enhance the performance of a computational algorithm, where nodes represent 
computations and edges represent communications, the nodes are divided equally among the processors so that the 
communications between them are minimized. 

The goal of a network partitioning algorithm is to divide any given network into approximately equal size groups 
irrespective of node similarities. Community detection on the other hand finds groups that either have an inherent or 
an externally specified notion of similarity among nodes within groups. Furthermore, the number of communities in 
a network and their sizes are not known beforehand and they are established by the community detection algorithm. 

Many algorithms have been proposed to find community structures in networks. Hierarchical methods divide 
networks into communities, successively,based on a dis-similarity measure, leading to a series of partitions from the 
entire network to singleton communities [H EB| • Similarly one can also successively group together smaller communities 
based on a similarity measure leading again to a series of partitions [IH [23[ . Due to the wide range of partitions, 
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structural indices that measure the strength of community structures are used in determining the most relevant ones. 
Simulation based methods are also often used to find partitions with a strong community structure 10, 24], Spectral 
[l7L I25I ] and flow maximization (cut minimization) methods 0, H(| have been successfully used in dividing networks 
into two or more communities. 

In this paper, we propose a localized community detection algorithm based on label propagation. Each node is 
initialized with a unique label and at every iteration of the algorithm, each node adopts a label that a maximum 
number of its neighbors have, with ties broken uniformly randomly. As the labels propagate through the network 
in this manner, densely connected groups of nodes form a consensus on their labels. At the end of the algorithm, 
nodes having the same labels are grouped together as communities. As we will show, the advantage of this algorithm 
over the other methods is its simplicity and time efficiency. The algorithm uses the network structure to guide its 
progress and does not optimize any specific chosen measure of community strengths. Furthermore, the number of 
communities and their sizes are not known a priori and are determined at the end of the algorithm. We will show that 
the community structures obtained by applying the algorithm on previously considered networks, such as Zachary's 
karate club friendship network and the US college football network, are in agreement with the actual communities 
present in these networks. 



II. DEFINITIONS AND PREVIOUS WORK 



As mentioned earlier, there is no unique definition of a community. One of the simplest definitions of a community 
is a clique, that is, a group of nodes where there is an edge between every pair of nodes. Cliques capture the intuitive 
notion of a community [6| where every node is related to every other node and hence have strong similarities with each 
other. An extension of this definition was used by Palla et al in [l4j, who define a community as a chain of adjacent 
cliques. They define two fc-cliques (cliques on fc nodes) to be adjacent if they share fc — 1 nodes. These definitions are 
strict in the sense that the absence of even one edge implies that a clique (and hence the community) no longer exists, 
fc-clans and fc-clubs are a more relaxed definitions while still maintaining a high density of edges within communities 
[3]. A group of nodes is said to form a fc-clan if the shortest path length between any pair of nodes, or the diameter 
of the group, is at most fc. Here the shortest path only uses the nodes within the group. A fc-club is defined similarly, 
except that the subnetwork induced by the group of nodes is a maximal subgraph of diameter k in the network. 

Definitions based on degrees (number of edges) of nodes within the group relative to their degrees outside the group 
were given by Radicchi et al If df and d° ut are the degrees of node i within and outside of its group U, then U 
is said to form a strong community if df 1 > d° ut , Vz G U. If YlieU^i 1 > J^ icrr d° ut , then U is a community in the 
weak sense. Other definitions based on degrees of nodes can be found in [6(. 

There can exist many different partitions of nodes in the network that satisfy a given definition of community. In 
most cases [3, [H, HI, E3, HI], the groups of nodes found by a community detection algorithm are assumed to be 
communities irrespective of whether they satisfy a specific definition or not. To find the best community structures 
among them we need a measure that can quantify the strength of a community obtained. One of the ways to measure 
the strength of a community is by comparing the density of edges observed within the community with the density 
of edges in the network as a whole [f|. If the number of edges observed within a community U is ejj, then under 
the assumption that the edges in the network are uniformly distributed among pairs of nodes, we can calculate the 
probability P that the expected number of edges within U is larger than eu- If P is small, then the observed density in 
the community is greater than the expected value. A similar definition was recently adopted by Newman [l3|, where 
the comparison is between the observed density of edges within communities and the expected density of edges within 
the same communities in randomized networks that nevertheless maintain every node's degree. This was termed the 
modularity measure Q, where Q = J2i( e a ~ a f), e a ^ the observed fraction of edges within group i and a\ is the 
expected fraction of edges within the same group i. Note that if is the fraction of edges in the network that run 
between group i and group j, then = J2j e ij- Q = implies that the density of edges within groups in a given 
partition is no more than what would be expected by a random chance. Q closer to 1 indicates stronger community 
structures. 

Given a network with n nodes and m edges N(n,m), any community detection algorithm finds subgroups of 
nodes. Let C\, C'2, C p be the communities found. In most algorithms, the communities found satisfy the following 
constraints 

1. C t nCj ■ = for i and 

2. [J. d spans the node set in N 

A notable exception is Palla et al [l4| who define communities as a chain of adjacent fc-cliques and allow community 
overlaps. It takes exponential time to find all such communities in the network. They use these sets to study 
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FIG. 1: An illustration of a dendrogram which is a tree representation of the order in which nodes are segregated into different 
groups or communities. 



the overlapping structure of communities in social and biological networks. By forming another network where a 
community is represented by a node and edges between nodes indicate the presence of overlap, they show that such 
networks are also heterogeneous (fat-tailed) in their node degree distributions. Furthermore, if a community has 
overlapping regions with two other communities, then the neighboring communities are also highly likely to overlap. 

The number of different partitions of a network N(n, m) into just two disjoint subsets is 2 n and increases exponen- 
tially with n. Hence we need a quick way to find only relevant partitions. Girvan and Newman Q proposed a divisive 
algorithm based on the concept of edge betweenness centrality , that is, the number of shortest paths among all pairs 
of nodes in the network passing through that edge. The main idea here is that edges that run between communities 
have high betweenness values than those that lie within communities. By successively recalculating and removing 
edges with highest betweenness values, the network breaks down into disjoint connected components. The algorithm 
continues until all edges are removed from the network. Each step of the algorithm takes 0(mn) time and since there 
are m edges to be removed, the worst case running time is 0(m 2 n). As the algorithm proceeds one can construct a 
dendrogram (see figure [I} depicting the breaking down of the network into disjoint connected components. Hence for 
any given h such that 1 < h < n, at most one partition of the network into h disjoint subgroups is found. All such 
partitions in the dendrogram are depicted, irrespective of whether or not the subgroups in each partition represent a 
community. Radicchi et al [la ] propose another divisive algorithm where the dendrograms are modified to reflect only 
those groups that satisfy a specific definition of a community. Further, instead of edge betweenness centrality, they 
use a local measure called edge clustering coefficient as a criterion for removing edges. The edge clustering coefficient 
is defined as the fraction of number of triangles a given edge participates in, to the total number of possible such 
triangles. The clustering coefficient of an edge is expected to be the least for those running between communities 
and hence the algorithm proceeds by removing edges with low clustering coefficients. The total running time of this 
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divisive algorithm is O(^). 

Similarly one can also define a topological similarity between nodes and perform an agglomerative hierarchical 
clustering [23l [2^|. In this case, we begin with nodes in n different communities and group together communities 
that are the most similar. Newman [12] proposed an amalgamation method (similar to agglomerative methods) using 
the modularity measure Q, where at each step those two communities are grouped together that give rise to the 
maximum increase or smallest decrease in Q. This process can also be represented as a dendrogram and one can cut 
across the dendrogram to find the partition corresponding to the maximum value of Q (see figure [T]) . At each step of 
the algorithm one compares at most m pairs of groups and requires at most 0(n) time to update the Q value. The 
algorithm continues until all the n nodes are in one group and hence the worst case running time of the algorithm 
is 0(n(m + nj). The algorithm of Clauset et al [3(| is an adaptation of this agglomerative hierarchical method, but 
uses a clever data structure to store and retrieve information required to update Q. In effect, they reduce the time 
complexity of the algorithm to 0(md log n), where d is the depth of the dendrogram obtained. In networks that 
have a hierarchical structure with communities at many scales, d ~ log n. There have also been other heuristic and 
simulation based methods that find partitions of a given network maximizing the modularity measure Q (ToL HH ■ 

Label flooding algorithms have also been used in detecting communities in networks [Ull^]. In (27[, the authors 
propose a local community detection method where a node is initialized with a label which then propagates step by 
step via the neighbors until it reaches the end of the community, where the number of edges proceeding outward from 
the community drops below a threshold value. After finding the local communities at all nodes in the network, an 
n X n matrix is formed, where the ij th entry is 1 if node j belongs to the community started from i and otherwise. 
The rows of the matrix are then rearranged such that the similar ones are closer to each other. Then, starting from 
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the first row they successively include all the rows into a community until the distance between two successive rows 
is large and above a threshold value. After this a new community is formed and the process is continued. Forming 
the rows of the matrix and rearranging them requires 0(n 3 ) time and hence the algorithm is time-consuming. 

Wu and Huberman (26j propose a linear time (0(m + n)) algorithm that can divide a given network into two 
communities. Suppose that one can find two nodes (x and y) that belong to two different communities, then they 
are initialized with values 1 and respectively. All other nodes are initialized with value 0. Then at each step of 
the algorithm, all nodes (except x and y) update their values as follows. If Zi, Z2, Zk are neighbors of a node z, 
then the value V z is updated as Vz i +Vz 2+---+ V 'k _ This process continues until convergence. The authors show that the 
iterative procedure converges to a unique value, and the convergence of the algorithm does not depend on the size n 
of the network. Once the required convergence is obtained, the values are sorted between and 1. Going through 
the spectrum of values in descending order, there will be a sudden drop at the border of two communities. This gap 
is used in identifying the two communities in the network. A similar approach was used by Flake et al |9| to find 
the communities in the WWW network. Here, given a small set of nodes (source nodes), they form a network of 
webpages that are within a bounded distance from the sources. Then by designating (or artificially introducing) sink 
nodes, they solve for the maximum flow from the sources to the sinks. In doing so one can then find the minimum 
cut corresponding to the maximum flow. The connected component of the network containing the source nodes after 
the removal of the cut set is then the required community. 

Spectral bisection methods (2o| have been used extensively to divide a network into two groups so that the number 
of edges between groups is minimized. Eigenvectors of the Laplacian matrix (L) of a given network are used in the 
bisection process. It can be shown that L has only real non-negative eigenvalues (0 < Ai < A2 < ... < A n ) and 
minimizing the number of edges between groups is same as minimizing the positive linear combination M = sf Ai, 
where Si — uf z and Uj is the eigenvector of L corresponding to A^. z is the decision vector whose z th entry can 
be either 1 or -1 denoting which of the two groups node i belongs to. To minimize M, z is chosen as parallel as 
possible to the eigenvector corresponding to the second smallest eigenvalue (The smallest eigenvalue is and choosing 
z parallel to the corresponding eigenvector gives a trivial solution). This bisection method has been extended to 
finding communities in networks that maximize the modularity measure Q [25| . Q can be written as a positive linear 
combination of eigenvalues of the matrix B, where B is defined as difference of the two matrices A and P. Ay is the 
observed number of edges between nodes i and j and Pij is the expected number of edges between i and j if the edges 
fall randomly between nodes, while maintaining the degree of each node. Since Q has to be maximized, z is chosen 
as parallel as possible to the eigenvector corresponding to the largest eigenvalue. 

Since many real-world complex networks are large in size, time efficiency of the community detection algorithm is an 
important consideration. When no a priori information is available about the likely communities in a given network, 
finding partitions that optimize a chosen measure of community strength is normally used. Our goal in this paper is 
to develop a simple time-efficient algorithm that requires no prior information (such as number, sizes or central nodes 
of the communities) and uses only the network structure to guide the community detection. The proposed mechanism 
for such an algorithm which does not optimize any specific measure or function is detailed in the following section. 

III. COMMUNITY DETECTION USING LABEL PROPAGATION 

The main idea behind our label propagation algorithm is the following. Suppose that a node x has neighbors 
xi, X2, Xk and that each neighbor carries a label denoting the community to which they belong to. Then x 
determines its community based on the labels of its neighbors. We assume that each node in the network chooses to 
join the community to which the maximum number of its neighbors belong to, with ties broken uniformly randomly. 
We initialize every node with unique labels and let the labels propagate through the network. As the labels propagate, 
densely connected groups of nodes quickly reach a consensus on a unique label (see figure [2]). When many such dense 
(consensus) groups are created throughout the network, they continue to expand outwards until it is possible to do 
so. At the end of the propagation process, nodes having the same labels are grouped together as one community. 

We perform this process iteratively, where at every step, each node updates its label based on the labels of its 
neighbors. The updating process can either be synchronous or asynchronous. In synchronous updating, node x at 
the t th iteration updates its label based on the labels of its neighbors at iteration t — 1. Hence, C x (t) = f(C Xl (t — 
1), C Xk (t — 1)), where c x (t) is the label of node x at time t. The problem however is that subgraphs in the 
network that are bi-partite or nearly bi-partite in structure lead to oscillations of labels (see figure [3]). This is 
especially true in cases where communities take the form of a star graph. Hence we use asynchronous updating where 
C x {t) = f{C Xil (*)>•••; C Xim (t), C Xi(m+1) (t — 1), C Xik (t — 1)) and Xn, Xi m are neighbors of x that have already been 
updated in the current iteration while Xj( m+ n, ...,Xik are neighbors that are not yet updated in the current iteration. 
The order in which all the n nodes in the network are updated at each iteration is chosen randomly. Note that while 
we have n different labels at the beginning of the algorithm, the number of labels reduces over iterations, resulting in 
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FIG. 2: Nodes are updated one by one as we move from left to right. Due to a high density of edges (highest possible in this 
case), all nodes acquire the same label. 
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FIG. 3: An example of a bi-partite network in which the label sets of the two parts are disjoint. In this case, due to the choices 
made by the nodes at step t, the labels on the nodes oscillate between a and b. 



only as many unique labels as there are communities. 

Ideally the iterative process should continue until no node in the network changes its label. However, there could be 
nodes in the network that have an equal maximum number of neighbors in two or more communities. Since we break 
ties randomly among the possible candidates, the labels on such nodes could change over iterations even if the labels 
of their neighbors remain constant. Hence we perform the iterative process until every node in the network has a 
label to which the maximum number of its neighbors belong to. By doing so we obtain a partition of the network into 
disjoint communities, where every node has at least as many neighbors within its community as it has with any other 

Q . 

community. If C\, C p are the labels that are currently active in the network and 3 is the number of neighbors 
node i has with nodes of label Cj, then the algorithm is stopped when for every node i, 

If i has label C m then df m > df J Vj 

At the end of the iterative process nodes with the same label are grouped together as communities. Our stop 
criterion characterizing the obtained communities is similar (but not identical) to the definition of strong communities 
proposed by Radicchi et al [l5j . While strong communities require each node to have strictly more neighbors within 
its community than outside, the communities obtained by the label propagation process require each node to have 
at least as many neighbors within its community as it has with each of the other communities. We can describe our 
proposed label propagation algorithm in the following steps. 

1. Initialize the labels at all nodes in the network. For a given node x, C X (Q) = x. 

2. Set t = 1. 

3. Arrange the nodes in the network in a random order and set it to X . 

4. For each x G X chosen in that specific order, let C x (t) — f(C Xil (t), ■■■,C Xim (t),C x . (m+1) (t — 1), ...,C Xik (t — 1)). 
/ here returns the label occurring with the highest frequency among neighbors and ties are broken uniformly 
randomly. 

5. If every node has a label that the maximum number of their neighbors have, then stop the algorithm. Else, set 
t = t + 1 and go to (3). 

Since we begin the algorithm with each node carrying a unique label, the first few iterations result in various small 
pockets (dense regions) of nodes forming a consensus (acquiring the same label). These consensus groups then gain 
momentum and try to acquire more nodes to strengthen the group. However, when a consensus group reaches the 
border of another consensus group, they start to compete for members. The within-group interactions of the nodes can 
counteract the pressures from outside if there are less between-group edges than within-group edges. The algorithm 
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converges, and the final communities are identified, when a global consensus among groups is reached. Note that 
even though the network as one single community satisfies the stop criterion, this process of group formation and 
competition discourages all nodes from acquiring the same label in case of heterogeneous networks with an underlying 
community structure. In case of homogeneous networks such as Erdos - Renyi random graphs [3l| that do not have 
community structures, the label propagation algorithm identifies the giant connected component of these graphs as a 
single community. 

Our stop criterion is only a condition and not a measure that is being maximized or minimized. Consequently there 
is no unique solution and more than one distinct partition of a network into groups satisfies the stop criterion (see 
figures U and [5]). Since the algorithm breaks ties uniformly randomly, early on in the iterative process when possibilities 
of ties are high, a node may vote in favor of a randomly chosen community. As a result, multiple community structures 
are reachable from the same initial condition. 

If we know the set of nodes in the network that are likely to act as centers of attraction for their respective 
communities, then it would be sufficient to initialize such nodes with unique labels, leaving the remaining nodes 
unlabeled. In this case when we apply the proposed algorithm the unlabeled nodes will have a tendency to acquire 
labels from their closest attractor and join that community. Also, restricting the set of nodes initialized with labels 
will reduce the range of possible solutions that the algorithm can produce. Since it is generally difficult to identify 
nodes that are central to a community before identifying the community itself, here we give all nodes equal importance 
at the beginning of the algorithm and provide them each with unique labels. 
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FIG. 5: The grouping of US college football teams into conferences are shown in (a) and (i>). Each solution ((a) and (6)) is an 
aggregate of five different solutions obtained by applying the algorithm on the college football network. 



We apply our algorithm to the following networks. The first one is Zachary's karate club network which is a network 
of friendship among 34 members of a karate club 32]. Over a period of time the club split into two factions due to 
leadership issues and each member joined one of the two factions. The second network that we consider is the US 
college football network that consists of 115 college teams represented as nodes and has edges between teams that 
played each other during the regular season in the year 2000 Q . The teams are divided into conferences (communities) 
and each team plays more games within its own conference than inter-conference games. Next is the co-authorship 
network of 16726 scientists who have posted preprints on the condensed matter archive at www.arxiv.org; the edges 
connect scientists who co-authored a paper [331 ] - It has been shown that communities in co-authorship networks are 
made up by researchers working in the same field or are research groups 22]. Along similar lines one can expect 
an actor collaboration network to have communities containing actors of a similar genre. Here we consider an actor 
collaboration network of 374511 nodes and edges running between actors who have acted in at least one movie together 
Q. We also consider a protein-protein interaction network [34| consisting of 2115 nodes. The communities are likely 
to reflect functional groupings of this network. And finally we consider a subset of the world wide web (WWW) 
consisting of 325729 webpages within the nd.edu domain and hyperlinks interconnecting them Q. Communities here 
are expected to be groups of pages on similar topics. 



A. Multiple community structures 



Figure 0] shows three different solutions obtained for the Zachary's karate club network and figure [5] shows two 
different solutions obtained for the US college football network. We will show that even though we obtain different 
solutions (community structure), they are similar to each other. To find the percentage of nodes classified in the same 
group in two different solutions, we form a matrix M, where My is the number of nodes common to community i in one 
solution and community j in the other solution. Then we calculate f sam e = ^QZi maxjjMy} + £3 . max i{-&%})^p . 
Given a network whose communities are already known, a community detection algorithm is commonly evaluated 
based on the percentage (or number) of nodes that are grouped into the correct communities [22|, [26J. / same is similar, 
where by fixing one solution we evaluate how close the other solution is to the fixed one and vice versa. While f sa me 
can identify how close one solution is to another, it is however not sensitive to the seriousness of errors. For example, 
when few nodes from several different communities in one solution are fused together as a single community in another 
solution, the value of f S ame does not change much. Hence we also use Jaccard's index which has been shown to be 
more sensitive to such differences between solutions |35j . If a stands for the pairs of nodes that are classified in the 
same community in both solutions, b for pairs of nodes that are in the same community in the first solution and 
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FIG. 6: Similarities between five different solutions obtained for each network is tabulated. An entry in the i th row and j th 
column in the lower triangle of each of the table is the Jaccard's similarity index for solutions i and j of the corresponding 
network. Entries in the i th row and j column in the upper-triangle of the tables are the values of the measure fsame for 
solutions i and j in the respective networks. The range of modularity values Q obtained for the five different solutions is also 
given for each network. 



different in the second and c vice- versa, then Jaccard's index is defined as .? , . It takes values between and 1, 

' a+b+c 1 

with higher values indicating stronger similarity between the two solutions. Figure [B] shows the similarities between 
solutions obtained from applying the algorithm five different times on the same network. For a given network, the 
2 jth en ^ r y m lower triangle of the table is the Jaccard index for solutions i and j, while the ij th entry in the upper 
triangle is the measure f sa me for solutions i and j. We can see that the solutions obtained from the five different runs 
are similar, implying that the proposed label propagation algorithm can effectively identify the community structure 
of any given network. Moreover, the tight range and high values of the modularity measure Q obtained for the five 
solutions (figure [6]) suggest that the partitions denote significant community structures. 



B. Aggregate 



It is difficult to pick one solution as the best among several different ones. Furthermore, one solution may be able 
to identify a community that was not discovered in the other and vice-versa. Hence an aggregate of all the different 
solutions can provide a community structure containing the most useful information. In our case a solution is a set 
of labels on the nodes in the network and all nodes having the same label form a community. Given two different 
solutions, we combine them as follows; let C 1 denote the labels on the nodes in solution 1 and C 2 denote the labels 
on the nodes in solution 2. Then, for a given node x, we define a new label as C x = (C*, C 2 ) (see figure [7]). Starting 
with a network initialized with labels C we perform the iterative process of label propagation until every node in the 
network is in a community to which the maximum number of its neighbors belong to. As and when new solutions 
are available they are combined one by one with the aggregate solution to form a new aggregate solution. Note that 
when we aggregate two solutions, if a community T in one solution is broken into two (or more) different communities 
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FIG. 7: An example of aggregating two community structure solutions. ti,t2,ts and <4 are labels on the nodes in a network 
obtained from solution 1 and denoted as C 1 . The network is partitioned into groups of nodes having the same labels. si,S2 
and S3 are labels on the nodes in the same network obtained from solution 2 and denoted as C 2 . All nodes that had label ii 
in solution 1 are split into two groups with each group having labels si and S2 respectively. While all nodes with labels t-j or 
<4 or ts in solution 1 have labels S3 in solution 2. C represents the new labels defined from C 1 and C 2 . 



S\ and S2 in the other, then by defining the new labels as described above we are showing preferences to the smaller 
communities Si and S2 over T. This is only one of the many ways in which different solutions can be aggregated. For 
other methods of aggregation used in community detection refer to [26|, [36|, [37| ■ 

Figure [5] shows the similarities between aggregate solutions. The algorithm was applied on each network 30 times 
and the solutions were recorded. An ij th entry is the Jaccard index for the aggregate of the first hi solutions with the 
aggregate of the first 5j solutions. We observe that the aggregate solutions are very similar in nature and hence a 
small set of solutions (5 in this case) can offer as much insight about the community structure of a network as can a 
larger solution set. In particular, the WWW network which had low similarities between individual solutions (Jaccard 
index range 0.4883 - 0.5931), shows considerably improved similarities (Jaccard index range 0.6604 - 0.7196) between 
aggregate solutions. 

IV. VALIDATION OF THE COMMUNITY DETECTION ALGORITHM 

Since we know the communities present in Zachary's karate club and the US football network, we explicitly verify 
the accuracy of the algorithm by applying it on these networks. We find that the algorithm can effectively unearth 
the underlying community structures in the respective networks. The community structures obtained by using our 
algorithm on Zachary's karate club network is shown in figure [4] While all the three solutions are outcomes of the 
algorithm applied to the network, figure @J6) reflects the true solution [33] ■ 

Figure \5\ gives two solutions for the US college football network. The algorithm was applied to this network 10 
different times and the two solutions are the aggregate of the first five and remaining five solutions. In both figures[5](a) 
and[5j6), we can see that the algorithm can effectively identify all the conferences with the exception of Sunbelt. The 
reason for the discrepancy is the following: among the seven teams in the Sunbelt conference, four teams (Sunbelt4 
= {North- Texas, Arkansas State, Idaho, New Mexico State}) have all played each other and three teams (Sunbelt3 = 
{Louisiana-Monroe, Middle- Tennessee State, Louisiana-Lafayette}) have again played one another. There is only one 
game connecting Sunbelt4 and Sunbelt3, namely the game between North- Texas and Louisiana-Lafayette. However, 
four teams from the Sunbelt conference (two each from Sunbelt4 and S unbelts) have together played with seven 
different teams in the Southeastern conference. Hence we have the Sunbelt conference grouped together with the 
Southeastern conference in figure [5][a). In figure [5^6), the Sunbelt conference breaks into two, with Sunbelt3 grouped 
together with Southeastern and Sunbelt4 grouped with an independent team (Utah State), a team from Western 
Atlantic (Boise State), and the Mountain West conference. The latter grouping is due to the fact that every member 
of Sunbelt4 has played with Utah State and with Boise State, who have together played five games with four different 
teams in Mountain West. There are also five independent teams which do not belong to any specific conference and 
are hence assigned by the algorithm to a conference where they have played the maximum number of their games. 

V. TIME COMPLEXITY 

It takes a near-linear time for the algorithm to run to its completion. Initializing every node with unique labels 
requires 0(n) time. Each iteration of the label propagation algorithm takes linear time in the number of edges (0(m)). 
At each node x, we first group the neighbors according to their labels (0(d x )). We then pick the group of maximum 
size and assign its label to x, requiring a worst-case time of 0{d x ). This process is repeated at all nodes and hence 
an overall time is 0(m) for each iteration. 
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FIG. 8: Similarities between aggregate solutions obtained for each network. An entry in the i row and j 1 column in the tables 
is Jaccard's similarity index between the aggregate of the first 5i and the first 5j solutions. While similarities between solutions 
for the karate club friendship network and the protein-protein interaction network are represented in the lower triangles of the 
first two tables, the entries in the upper triangle of these two tables are for the US college football network and the co-authorship 
network respectively. The similarities between aggregate solutions for the WWW is given in the lower triangle of the third 
table. 



As the number of iterations increases, the number of nodes that are classified correctly increases. Here we assume 
that a node is classified correctly if it has a label that the maximum number of its neighbors have. From our 
experiments, we found that irrespective of n, 95% of the nodes or more are classified correctly by the end of iteration 
5. Even in the case of Erdos - Renyi random graphs [3l[ with n between 100 and 10000 and average degree 4, which 
do not have community structures, by iteration 5, 95% of the nodes or more are classified correctly. In this case, the 
algorithm identified all nodes in the giant connected component as belonging to one community. 

When the algorithm terminates it is possible that two or more disconnected groups of nodes have the same label (the 
groups are connected in the network via other nodes of different labels). This happens when two or more neighbors 
of a node receive its label and pass the labels in different directions, which ultimately leads to different communities 
adopting the same label. In such cases, after the algorithm terminates one can run a simple breadth-first search on 
the sub-networks of each individual groups to separate the disconnected communities. This requires an overall time 
of 0{m + n). When aggregating solutions however, we rarely find disconnected groups within communities. 



VI. DISCUSSION AND CONCLUSIONS 



The proposed label propagation process uses only the network structure to guide its progress and requires no external 
parameter settings. Each node makes its own decision regarding the community to which it belongs to based on the 
communities of its immediate neighbors. These localized decisions lead to the emergence of community structures in a 
given network. We verified the accuracy of community structures found by the algorithm using Zachary's karate club 
and the US college football networks. Furthermore, the modularity measure Q was significant for all the solutions 
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FIG. 9: The cumulative probability distributions of community sizes (s) are shown for the WWW, co-authorship and actor 
collaboration networks. They approximately follow power-laws with the exponents as shown. 



obtained, indicating the effectiveness of the algorithm. Each iteration takes a linear time O(m), and although one 
can observe the algorithm beginning to converge significantly after about 5 iterations, the mathematical convergence 
is hard to prove. Other algorithms that run in a similar time-scale include the algorithm of Wu and Huberman [26| 
(with time complexity 0(m + n j) and that of Clauset et al [3(| which has a running time of (9(nlog 2 n). 

The algorithm of Wu and Huberman is used to break a given network into only two communities. In this iterative 
process two chosen nodes are initialized with scalar values 1 and and every node updates its value as the average of 
the values of its neighbors. At convergence, if a maximum number of a node's neighbors have values above a given 
threshold then so will the node. Hence a node tends to be classified to a community to which the maximum number of 
its neighbors belong. Similarly if in our algorithm we choose the same two nodes and provide them with two distinct 
labels (leaving the others unlabeled), the label propagation process will yield similar communities as the Wu and 
Huberman algorithm. However to find more than two communities in the network, the Wu and Huberman algorithm 
needs to know a priori how many communities there are in the network. Furthermore, if one knows that there are 
c communities in the network, the algorithm proposed by Wu and Huberman can only find communities that are 
approximately of the same size, that is — , and it is not possible to find communities with heterogeneous sizes. The 
main advantage of our proposed label propagation algorithm over the Wu and Huberman algorithm is that we do not 
need a priori information on the number and sizes of the communities in a given network; indeed such information 
usually is not available for real-world networks. Also, our algorithm does not make restrictions on the community 
sizes. It determines such information about the communities by using the network structure alone. 

In our test networks, the label propagation algorithm found communities whose sizes follow approximately a power- 
law distribution P(S > s) ~ s~ u with the exponent v ranging between 0.5 and 2 (figure [9|) . This implies that there 
is no characteristic community size in the networks and it is consistent with previous observations [22l . [30L l38j . While 
the community size distributions for the WWW and co-authorship networks approximately follow power-laws with a 
cut-off, with exponents 1.15 and 1.98 respectively, there is a clear crossover from one scaling relation to another for 
the actor collaboration network. The community size distribution for the actor collaboration network has a power-law 
exponent of 2 for sizes up to 164 nodes and 0.5 between 164 and 7425 nodes (see figure [S])- 

In the hierarchical agglomerative algorithm of Clauset et al [30l | , the partition that corresponds to the maximum Q 
is taken to be the most indicative of the community structure in the network. Other partitions with high Q values 
will have a structure similar to that of the maximum Q partition, as these solutions are obtained by progressively 
aggregating two groups at a time. Our proposed label propagation algorithm on the other hand finds multiple 
significantly modular solutions that have some amount of dissimilarity. For the WWW network in particular, the 
similarity between five different solutions is low, with the Jaccard index ranging between 0.4883 to 0.5921, yet all 
five are significantly modular with Q between 0.857 to 0.864. This implies that the proposed algorithm can find not 
just one but multiple significant community structures, supporting the existence of overlapping communities in many 
real- world networks 1141. 
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